Finite temperature fidelity susceptibility for one-dimensional quantum systems 
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Wo calculate the fidelity susceptibility Xf for the Luttingcr model and show that there is a uni- 
versal contribution linear in temperature T (or inverse length l/L). Furthermore, we develop an 
algorithm - based on a lattice path integral approach - to calculate the fidelity F{T) in the ther- 
modynamic limit for one-dimensional quantum systems. We check the Luttinger model predictions 
by calculating Xfi^) analytically for free spinless fermions and numerically for the XX Z chain. 
Finally, we study Xf at the two phase transitions in this model. 
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Phase transitions are usually identified by considering 
suitably defined order parameters. Lately, new concepts 
originating from quantum information theory have been 
put forward which allow to detect phase transitions with- 
out any prior knowledge of the order parameter [1-12]. 
The most widely used measures are the entanglement en- 
tropy [1] and the fidelity [2-10]. The latter approach is 
based on the notion that at a quantum phase transition 
the ground state wave function is expected to change dra- 
matically with respect to a parameter A driving the tran- 
sition [5]. If the Hamiltonian is given by H\ = Ho + AO, 
then the fidelity is defined as 

Fo{X) = V(*o|*a)(*a|*o)/(*o|*o)(*a|*a) (1) 

where \^o) [I* a)] is the ground state wave function of Hq 
[H\], respectively. The fidelity has been studied analyti- 
cally for one-dimensional (ID) models like the transverse 
Ising or the XY model [5, 7, 8] as well as numerically 
for a number of other systems [2, 3, 8]. Importantly, 
the fidelity approach connects many different areas of 
physics and is not restricted to the study of phase tran- 
sitions. The overlap between wave functions also plays a 
central role for scattering problems (Anderson's orthog- 
onality catastrophy) [13], as a measure for variational 
wave functions, for quantum information processing [14], 
the Loschmidt echo [15], and for quench dynamics [16]. 
Apart from calculating the fidelity for specific models it 
is therefore of great interest to understand possible uni- 
versal behavior. For critical ID quantum systems such 
universality is often related to conformal invariance. Im- 
portant examples arc the scaling of the free energy [17] 
and the entanglement entropy [1] with system size L and 
temperature T. 

In this letter we will introduce a new finite temperature 
(mixed state) fidelity and show that it leads to the fidelity 
susceptibility Xf used in recent quantum Monte Carlo 
simulations [3]. We then show that Xf for the Luttingcr 
model has a universal term linear T. Similarly, there is 
a universal term ~ 1/i for a finite system at zero tem- 
perature. Xj{T = 0) in the thermodynamic limit, on the 
other hand, depends on a cutoff, a fact, which has been 



missed in an earlier work [10]. Furthermore, we express 
F{T) in the thermodynamic limit for any ID quantum 
system as a function of the largest eigenvalues of three 
transfer matrices. This allows for a very efficient numer- 
ical calculation of the fidelity making it an ideal tool for 
finding phase transitions without any prior knowledge of 
the order parameters. We apply this method to study 
X/(T) for the 5 = 1/2 XX Z chain with respect to a 
small change in the anisotropy A allowing us to check 
om results for the Luttinger model directly. A further 
check is provided by an analytic calculation of Xf{T) in 
the free fermion case. Finally, we extract Xf{T = 0) for 
the XX Z model from the numerical data and discuss its 
behavior at the two critical points. 

We can generalize (1) to finite temperatures so that 
Ft{0) = 1 and limT^o^rCA) = Fo{X) by 

Ft (A) = ^TV{e-^^o/2e-0HA/2}/(^o z^)i/4 (2) 

where /3 = 1/T, = Tre-'^^o, and Zx = Tre"''^^. For 
a many-body system the fidelity is expected to vanish 
exponentially with the number of particles N no matter 
how small the driving parameter A is [13]. The fidelity 
density /(A) = —^\nF, however, stays finite. Since 
/(A = 0) = is a minimum, the first term in an expan- 
sion for small A vanishes giving rise to the definition of 
the fidelity susceptibifity Xf = i.d'^ f / d\^)\=o [9]. Prom 
Eq. (2) we find that 

X/ = ^ j^'\dT{(J-d{r)dm - {Of} (3) 

where T denotes time ordering and 0{t) = 

exp(TiJo)0 cxp(— riJo)- In the following, we will 
consider the case 0(r) = X]^o(r, r) where 6(r, r) is a 
local operator. By using a Lehmann representation, 
Eq. (3) can be shown to be consistent for T — > with 
the ground state fidelity directly obtained from the 
definition (1) [9]. Eq. (3) has previously been used to 
define Xf{T) [3]. Here this expression for Xf{T) in terms 
of a correlation function directly follows from Eq. (2). 



2 



Note, however, that F{T) in (2) is different from the 
mixed state fideUty as defined in [6, 7] which does 
not allow to express the corresponding %/ as a simple 
correlation function. Importantly, it has been shown 
that if Xf as obtained from the mixed state fidelity 
in [6, 7] diverges then so does Xf as given in (3) and 
vice versa [3]. Finally, wc note that if [Hq,0] — then 
X/(T) = x/8T with X = ((E.o.)2)/(iVr) being the 
regular susceptility. 

The generic low-energy effective theory for a gapless 
ID quantum system is the Luttinger model [18] 



L/2 



dx 



L/2 



K 



(4) 



Here w is a velocity, L = Na the length with a being the 
lattice constant, and K the Luttinger parameter. is 
a bosonic field obeying the standard commutation rule 



(x),n(x')] = iSix-x') with n 



dr4'- In general. 



both K and v will change as a function of a driving pa- 
rameter A in the Hamiltonian of the microscopic model. 

The operator appearing in (3) is therefore given by 
O = 6i 62 with 



Ok 
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L/2 



dx 



L/2 



K 



K 



(9.0)' 



(5) 



and ai = dv/dX, 012 = v{dK / dX) / K . We note that Oi 
is proportional to the Hamiltonian itself. By rescaling 

we can express the Hamilto- 
nian and therefore also 0\ as the sum of the holo- and an- 
tiholomorphic components of the energy-momentum ten- 
sor [19]. The finite temperature correlation function (3) 
involving 0\ can then be calculated with the help of the 
operator product expansion for this conformally invari- 
ant theory. While the cross term vanishes, the integral 
(3) for the operator O2 is divergent and we introduce 
a cutoff by replacing j^^"^ It/"^' Combining both 
contributions we find in the thermodynamic limit at low 
temperatures 



with A = l/{'KVTo) and c = 1 being the central charge 
of the free bosonic model. The universality found here 
for the leading linear temperature dependence of Xf is 
reminiscent of the universal term in the free energy of 
ID critical quantum systems quadratic in temperature 
[17]. We also want to remark that a universal subleading 
term in the zero temperature fidelity has recently been 
discovered in certain systems [20]. 

X/(T = 0) as obtained in (6), on the other hand, is 
cutoff dependent. This seems to be in contrast to an 
earlier work [10] where Xf "^as directly calculated at 
zero temperature using the definition (1). This leads to 
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FIG. 1: Transfer matrices for calculating F{T). Each 
open [shaded] plaquette represents a local Boltzmann weight 
exp{—ehl^''~~^^) [oxp(— e/i^'''"'"'^)], respectively, with e being the 
Trotter parameter. 



Xf = {dK/dX)^/{4:NK^)J2k>o and the result in [10] is 
obtained if one assumes N/2 fc-valucs in the sum. The 
Luttinger model, however, is a continuum model and the 
sum therefore not restricted. If we introduce a UV cutoff 
NA/2 then the first term in (6) is reproduced. 

Similarly, we can calculate Xf for the Luttinger model 
of finite size L at zero temperature using Eq. (3). Due to 
the unusual imaginary-time integration the result cannot 
be obtained by simply replacing v/T ^ L but rather the 
second term in (6) gets replaced by c{dv/dX)^/{8v^L). 

By using a lattice path integral representation, a ID 
quantum model can be mapped onto a two-dimensional 
classical model with the additional dimension corre- 
sponding to the inverse temperature. For the fidelity 
(2) this amounts to separate Trotter-Suzuki decomposi- 
tions for each of the exponentials. We consider a Hamil- 
tonian with nearest-neighbor interaction and decompose 



the Hamiltonian into H^ 

Erodd^oA^^- "^^^^ allows US to write exp{—/3Ho) = 
limM->oo[exp(— ei?o) exp(— effg)]*^ and equivalently for 
the other exponentials in (2). Here e = /3/M is the Trot- 
ter parameter. Rearranging the local Boltzmann weights 
we can define the column transfer matrices depicted in 
Fig. 1. The spectra of these transfer matrices have a gap 
between the largest and the next-leading eigenvalue thus 
allowing it to perform the thermodynamic limit exactly 
[21]. For the fidelity density we find 



/-r(A) = -j-lnF = -^lnf ,X 



(7) 



where Ay, Ao, and A\ arc the largest eigenvalues of 
the transfer matrices T/, Tq, and Tx defined in Fig. 1, 
respectively. Because /t(0) = 9/t/5A|a=o = we 
can calculate the fidelity susceptibility by Xf{T) = 
2 limA->.o /t(A)/A^, i.e., without having to resort to nu- 
merical derivatives. The transfer matrices can be ef- 
ficiently extended in imaginary time direction - corre- 
sponding to a successive reduction in temperature - by 
using a density-matrix renormalization group algorithm 
applied to transfer matrices (TMRG). If we are mainly 
interested in Xf then only small parameters A have to be 
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FIG. 2: TMRG data (dashed) and fits (solid lines) for Xf{T) 
in the ferrmomagnetic Ising regime. Inset: TMRG data for 
the free fermion case (circles) compared to the exact solution 
(dashed line) and the low-T asymtotics (solid line). 



FIG. 3: TMRG data (dashed) and low temperature fits (solid 
lines) for Xf with the slope fixed by (6). The main fig- 
ure (inset) shows data for A — — 0.8, — 0.7, • • • ,—0.1 (A = 
0.1, 0.2, • • • , 0.6) from top to bottom, respectively. 



considered, allowing it to renormalize all three transfer 
matrices with the same reduced density matrix. Apart 
from the two different Boltzmann weights necessary to 
form the three transfer matrices depicted in Fig. 1 the al- 
gorithm can therefore proceed in exactly the same way as 
the TMRG algorithm to calculate thermodynamic quan- 
tities. For technical details of the algorithm the reader is 
therefore referred to Refs. [21, 22]. 

In the following, we want to study XfiT) for the XX Z 
model defined by 

r 

with respect to a change in anisotropy A. Here S is 
a spin 5' = 1/2 operator and J the exchange constant 
which we set to 1. We note that Xf zero temperature 
for finite chains has previously been studied in [2, 8]. The 
model is gap less for —1 < A < 1 and gapped for |A| > 1. 
At A = the model describes non-interacting spinless 
fermions and Xf can be calculated exactly. The various 
diagrams can be combined into two contributions 
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with X = cos ki + cos k2 
cos2(fci — fca) — cos(A:i - 



y 



cos fca 
«3j - 



cos(fci 



k. 



■k3)cos{k2-k3), nf 



- h), 

2-K3), = 1/[1 + 

exp(/3cosfc)] and = 1 — . The first contribution 
at low temperatures yields xf'' = 0.19537(±5) + ©(T^) 

whereas the second is given hy Xf — T/{6n). In the in- 
set of Fig. 2 the exact solution for A = is compared with 
the TMRG data obtained from Xf = 2/(A + 6A)/{6A)^ 
with SA = 10"^. The relative error without any extrap- 
olation is less than 0.1% for T > 0.1 and less than 1% for 
T > 0.04. 



In the gapped regime, |A| > 1, the fidelity suscepti- 
bility will show activated behavior. Following the argu- 
ments in [23] for the magnetic susceptibility we expect 
Xf T~3/2 exp(— 7/T) with 7 = — A — 1 being the spec- 
tral gap for A < — 1 and 7 being half the spectral gap for 
A > 1. Note that in the latter case the spectral gap is 
exponentially small for A > 1 making it difficult to de- 
tect numerically. As shown in Fig. 2 a fit of the data for 
A < — 1 is consistent with this scaling form with fitted 7 
values close to the one theoretically expected. 

In the gapless regime, — 1 < A < 1, we know from 
the Bethe ansatz that K = 7r/(7r — arccosA) and v = 
TT\/l — A2/(2 arccos A). This allows us to check the uni- 
versality of the contribution linear in T in (6) by a direct 
comparison with the TMRG data (see Fig. 3) and to ac- 
curately extract x/(^ = 0)- This method fails, however, 
close to A = 1 where corrections due to Umklapp scat- 
tering start to become important (as will be discussed 
in more detail below) as well as very close to A = — 1 
where the Luttinger model fails because the dispersion 
of the elementary excitations becomes quadratic. The 
fidelity susceptibility as a function of A for various tem- 
peratures as well as the extrapolated T = curve are 
shown in Fig. 4(a). Comparing with the theoretical re- 
sult (6) we can also extract the momentum cutoff A (see 
Fig. 4(b)). There is a clear divergence of Xf at the first 
order phase transition A = — 1. A fit of the extrapo- 
lated zero temperature curve shown in Fig. 4(a) gives 
X/(A > -1) - 0.017/(A + 1)1-26. This requires that 
the cutoff A vanishes for A —1 because otherwise we 
would find from (6) a divergence ~ 1/(A -I- 1)^ as pre- 
dicted in [10]. Indeed, a fit of the extracted momentum 
cutoff as shown in Fig. 4(b) yields A - 0.43(A + 1)°-^^ 
and therefore x/(A > -1) 0.013/(A + l)'^-^^ which is 
consistent with the direct fit. 

At the Kosterlitz-Thouless (KT) phase transition, A = 
1, on the other hand, the behavior is different. Here the 
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FIG. 4: (a) The solid line depicts the extrapolated T — 
curve in the critical regime, the symbols the finite temper- 
ature curves as indicated on the plot, (b) The momentum 
cutoff A as a function of A (symbols) . The dashed line shows 
a fit for A > -1. (c) Xf{T) for A = 1, and (d) for A = -1 
with the dashed line representing a low-T fit as discussed in 
the text. 



be reliably tested by comparing with the TMRG data. 
While the term in (10) should dominate at very low tem- 
peratures, subleading corrections might be of equal im- 
portance in the temperature range accessible numerically. 

To summarize, we have shown that the fidelity suscep- 
tibility for the Luttinger model has a universal term lin- 
ear in temperature or inverse length. Apart from being 
relevant for quantum phase transitions we believe that 
this result is also important to analyze sudden quantum 
quenches. Furthermore, we have introduced a numerical 
method to calculate the finite temperature fidelity in the 
thermodynamic limit for any ID quantum system with 
short range interactions. Finally, based on a RG treat- 
ment, we have predicted a In^ T divergence of Xf s-t the 
KT transition in the XX Z model. 

JS thanks I. Affleck and F. Alet for valuable discus- 
sions. This work was supported by the MATCOR school 
of excellence. 



finite temperature data show that a maximum in x/ at 
A > 1 exists which shifts to smaller A with decreasing 
temperature. The dependence of the cutoff A near A = 1 
seems to be consistent with A ^ Ai + (1 — A)" with a 
constant Ai and an exponent a both greater than zero. 
If this is indeed the case, we find from (6) that Xfi^^ 
l)^Ai/[167r2(l-A)]. 

Finally, we want to discuss the temperature depen- 
dence of Xf right at the phase transitions. For A — 
— 1, shown in Fig. 4(d), we find a divergence Xf ~ 
0.002(±l)r-2-5(±i) where the error is determined by a 
variation of the fit interval. As argued above, we also 
expect X/(T) to diverge for T — > and A = 1. The 
numerical data, shown in Fig. 4(c), however, do not eas- 
ily allow to extract the low temperature behavior. If 
we assume that limA-i-i A = Ai > then we can calcu- 
late the temperature dependence analytically as follows. 
At the isotropic point, Umklapp scattering is marginally 
irrelevant and the Luttinger parameter has to be re- 
placed by a running coupling constant i^T — > 1 + g\\(l)/2 
where I ~ InTo/T with a scale Tq of order J and 
3||(Z) = {2K*-2)/tiiTih[{2K*-2)l+atanh{{2K*-2)/gp] 
[18]. For T = To we have Z = and gj|(0) = ,gj| while for 
T -J> it follows that / oo and K = 1 + gi\/2 ^ K* 
where K* is the fix point value. For I large we can neglect 
the part cx gjj. We therefore obtain 



X/(A = 1,T) 



hm ; r 

A^i 32(1+511/2)2 

T<sTo 2Ai 2frr irr\ 



dg\\ 
dA 



(10) 



While this prediction resolves some confusion about the 
behavior of Xf ^-t a KT transition [2, 8-10] it cannot 
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